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ABSTRACT 

In this paper, we have investigated the geodesies of neutral particles near a five-dimensional 
charged black hole using a comparative approach. The effective potential method is used to 
determine the location of the horizons and to study radial and circular trajectories. This also 
helps us to analyze the stability of radial and circular orbits. The radius of the innermost stable 
circular orbits have also been determined. Contrary to the case of massive particles for which, the 
circular orbits may have up to eight possible values of specific radius, we find that the photons will 
only have two distinct values for the specific radii of circular trajectories. Finally we have used 
the dynamical systems analysis to determine the critical points and the nature of the trajectories 
for the timelike and null geodesies. 

Subject headings: Geodesic Motions; Charged black hole; Effective potential; Horizons; Stability; Dy- 
namical systems; Critical points 



1. Introduction 

With the development of string theory, the 
study of black holes in higher-dimensional space- 
time s have assumed fundamental importanc e (see 
e.g. lEmparan fc Reall I 12005 iKanti l 12004 . In 
these models, the minimal mass of a black hole 
can be much smaller than the Planck mass (10 19 
GeV). Such mini black holes could be produced 
in high energy experiments: in col liders like LHC 



within the classical solutions of vacuum Einstein 
equations, by an asymptotically flat black hole in 
higher dimensions, including the case of the rotat- 
ing black holes. The study of higher dimensional 
black holes have gained m omentum in the first 
decade of t his millennium dGibbons et al~ 2004 ; 



ISenl |2005t iFrolov et al~l l2003h . Static, spheri 



and in cosmic ray experiments ([Argvres et al 
19981 ). Estimations show that a TeV-size black 



cally symmetric exterior vacuum solutions of the 
braneworld mod els were first proposed by Dad 
hich and others (Da dhich et al. 2000; iMaartens 



2004), with the metric having the structure of 



hole in TeV gravity is small enough to be described 
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the Reissner-Nordstrom (RN) solution and inter- 
preted as a black hole with a tidal charge arising 
via gravitational effects from the fifth dimension. 
Static, spherically symmetric charged black hole 
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solutions in 5-dimensions are known to be param- 
eterized by their mass a nd electric charges and rep - 
resent the RN solution (jEmparan fc Reall 1 120081 ) 
Extensive studies in higher-dimensional space- 
times over the last few decades have led many au- 
thors to inve stigate the g e odesic motions in such 



spacetimes (IPage et al. I l2007t ICardoso et al. 



2006t iKonoplva et al. 1 120081 ) . The geodesic struc- 
ture of Schwarzschild AdS bl ack hole has bee n 
considered by Cruz and others ( Cruz et al. | |2005). 



Motion of massive particles around a rotating 
black hole in a branewor ld has been studied 
(|Abduiabbarov et al. llioToh and the principal null 
geodesies and circular photon orbits for the Sen 



black hole have been investigated (jHioki et al 
2008). The effective potentials for radial null 



geodesies in RN-dS and Kerr-d S spacetimes were 
analyzed (|Stuchlik et al. Ill99ll ) and bound orbits 
in a Kerr spacetime have been discussed by some 



authors ( Kraniotis 2004) . Compact calculation 



of the perihelion precession of Mercury has been 
done, based on the inversio n problem of hyperellit 



tic integrals by Kraniotis ([Kraniotis et al. 200 



Analytic solutions of the geodesic equations in 
various spacetime geometries in 4-dimcnsions, as 
well as a number o f cases in higher dimens ions 
has been obtained jHackmann et al. 2008al ibli3: 
Kagramanova et al. I [2006). In this work, we 
have investigated the radial and circular trajec- 
tories for photons and massive particles in a five- 
dimensional RN spacetime using a comparative 
approach and have determined the fixed points of 
the phase trajectories. For such a non-rotating 
charged black hole, the solutions are uniquely 
characterized by th eir mass, charge and the cos- 
mological constant (jGibbons et al 1 120021) . 

The paper is organized as follows: In Section 
II the line element and the horizon function has 
been defined and the horizons are located. The 
equations for the 5-dimensional geodesies in the 
region surrounding the black hole are determined 
in the next section. In Section IV, the trajecto- 
ries of the test particles are analyzed using the 
method of effective Newtonian potential approx- 
imation. This equation is used to study radial 
motion and the corresponding stability of radial 
trajectories. To study circular motion of particles 
and the stability of circular orbits, the equations 
have been transformed in terms of the inverse of 
the radial coordinate and the energy and angular 



momentum of massive particles have been deter- 
mined. The radius of the innermost stable circular 
orbit of massive particles and photons have been 
calculated. Contrary to the case of massive parti- 
cles, for which the orbits may have up to eight pos- 
sible values of the specific radius, it is found that 
photons can have circular trajectories for only two 
distinct values of specific radius of the orbits. In 
Section V, we have made an analysis of the time- 
like and null geodesies using the dynamical sys- 
tems approach. The fixed points are determined 
and the nature of the phase trajectories have been 
analyzed. The summary and the conclusions are 
presented in Section VI. 

2. Preliminaries 

Let us consider a 5-dimensional spacetime in 
the presence of a cosmological constant. The field 
equations are given by 



Gab — —^-(s)9ab + k^Tab 



(1) 



where qab is a 5-dimensional metric of signature 
(- + + + +), Gab is the 5-dimensional Einstein 
tensor, Tab represents the 5-dimensional energy- 
momentum tensor and A( 5 ) is the 5-dimensional 
cosmological constant. The spacetime is assumed 



to be one of constant curvature K 



For AdS 



geometry, K is negative and (3 = — 1, whereas for 
dS, K is positive so that /3 — 1. The radius of 
curvature 



/ = 



3/3 



A 



(5) 



(2) 



of the spacetime provides the length scale neces- 
sary to have a horizon. To simplify our notation 
we shall drop the subscript and henceforth denote 
A( 5 ) by A. The above equations are satisfied by 
the exterior metric of the black hole field, which 
has the form 



dS 2 



-f(r)dt* 



dr 2 



r 2 dQi 



(3) 



where, df2§ = {d6 2 + 



+ sin 2 (f>dtp 2 )) is 



the metric of the unit 3-sphere. For this static, 
spherically symmetric exterior vacuum solution of 
the Einstein equations in RN spacetime, the lapse 
function f(r) is defined as 
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Here, q and M are the charge and the mass of 
the black hole, respectively. For a given value of 
M, q and A, the horizon function A depends only 
on the radial coordinate r. In diagonal form, the 
exterior metric can be written as 



dS 2 



A 



-dt 2 



A 



-dr 2 



r 2 dn 2 3 . 



(5) 



The spacetime, has an intrinsic singularity at r = 
0, the nature of which depends on the choice of 
A and q. These singularities may either be black 
holes or naked singularities. In fact there is no 
guarantee that "black hole" candidates are in- 
deed black holes and some (or all) of " black hole" 
candidates could instead be naked singu larities 
(|Virbhadra 1 120091: IVirbhadra et al. I l2002h . Fur- 
ther, the net electric charge of astronomical ob- 
jects are constrained within a certain range if their 
exterior spacetimes are described by the Reissner- 
Nordstr om metric, w hich induces power-law po- 
tentials ( Iorio 1 120111 ). Here we assume that both 



A and q are chosen in such a way that the space- 
time do not h ave any spacelike naked singularity 
(Carroll 1 l2004h . The lapse function vanishes at the 
zeros of the equation 



for which 



Ar b - 6r 



A = 0, 



24ikfV - 



= 0. 



(6) 



(7) 



The horizons are located at the real, positive 
zeros of A, indicating the coordinate singularities. 
The effective potential also vanishes at the zeros of 
the horizon function A, indicating the location of 
these horizons. The variation of the effective po- 
tential V e ff with specific radius r/M in the case of 
radial motion of massive particles, in the field of 
the black hole for different range of values of r/M, 
are shown in Fig. Q] The location of the horizons 
can be easily identified from these figures. As ex- 
pected, we find that the field of a charged black 
hole in the de Sitter spacetime is characterized 
by the presence of all three horizons: the Cauchy 
horizon, the event horizon and the cosmological 
horizon, whereas there are only two horizons in 
the Anti de Sitter spacetime. 

Among the six distinct roots of equation © in 
the RN spacetime, only two are real for the AdS 
case, with a single one being positive, which gives 



us the radius of the event horizon of the black hole. 
This positive root is found as 



(242V3M 2 ) 122V3 



3? 



A!ft 



1/3 



+ 



3? 



32V3A 



(8) 



where 3i = 3 (2 A + 2VB 3 + A 2 ) 
with A = 8 - 24A/ 2 A + 3q 4 A 2 and B = 8M 2 A - 4. 

3. Five-dimensional Geodesies 



The geodesic equations are 

d 2 Z A (5) r A dzB dz ° 

d\ 2 BC d\ dX 



(9) 



where A is an affine parameter along the geodesic 
curve z A (X), ^Tg C are the 5-dimensional Christof- 
fel symbols of the second kind and z A are the 
coordinates of the 5-dimensional spacetime. The 
nature of trajectories of the test particles depend 
on their energies and angular momentum as well as 
the cosmological constant and the electric charge 
of the black hole. The mass of the black hole 
can be absorbed by rescaling the radial coordi- 
nate. We consider only neutral test particles, so 
that their traj ectories describe geod esies in the 
RN spacetime ( Chandrasekhar 1983). 

On account of the spherical symmetry, the mo- 
tion of the test particles can be determined by an- 
alyzing their motion on the equatorial plane, for 
which 6,<j) = tt/2. Substituting the explicit form 
of the bulk metric, we obtain the following set of 
geodesic equations 



d 2 t B(r) dt dr 
dA 2 + ~Mr) dXdX 



= 0, 



d 2 r ., , , / dt 
— +A(r)B( r )(- 



B(r) 
A{r) 



(10) 



0, 



d 2 ip 1 dr dip 
~d>? + rdXdX 
where we have defined A(r) = 



0. 



(11) 
(12) 



-f(r) and B(r) 



Arf 
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We now proceed to make a study of the 
geodesies using a comparative approach. The 
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effective Newtonian potentiai formalism is used 
to determine the location of the horizons and to 
study radial and circular trajec tories. Most of the 



pap ers mentioned earl i er from (|Cruz et al. 1 12005) 
to ( Hackmann et al. I l2008a ) use this method. 



Subsequently, we determine the critical points and 
the nature of the trajectories of the particles us 



ing the dynamical system s analysis (jUzan et al 



2001; Dah ia et al. I 2007) 



4. Analysis of particle trajectories using 
the effective Newtonian potential ap- 
proach 

The lagrangian for particle motion in this RN- 
AdSs is given by 



where, e = 1 for timelike geodesies and e = for 
null geodesies, we have for the geodesies on the 
equatorial plane 



r 4 V e +7T J =^+^(r)( e - 



V — 

(18) 

Here we have used the sig n convention o f Carroll 
(jCarroll Il200ll) and Hartle (|Hartle Il2003h . We can 
rewrite (IT8t as 



1 f dr 



2 V dX 



E eff ~ V ef f(r), 



(19) 



where the equatorial geodesies defined above are 
subjected to an effective potential V e ff(r) given 

by 



g AB z A z B 



+ - 



Ar 2 
~6~ 



1 



Ar 2 
6 



2q12 + sin 2 d sin 2 ^ 2 ), (13) 



r 2 (6 2 + sin 



where, an overdot represents differentiation with 
respect to the affine parameter A along the 
geodesic z A (X). Since t and ip are cyclic coordi- 
nates, there are two conserved quantities, namely 
the energy E and the momentum conjugate to ip. 
Thus, we have 



„ dt „ . . dt d , , dt 



and 



dl 2 , 

p^p = — r = 2r ip = constant = 2L. 

dip 



so that 



L = r 



<bp_ 
dX' 



(14) 



(15) 



(16) 



where L is now the total angular momentum of the 
particles confined to the equatorial plane. From 
the normalization condition 



9ab- 



dz A dz B 
dX dX 



(17) 



Ve ff (r) 



A 



e + 




and 



Thus (fTTH) is the equation of motion of a parti- 
cle of unit mass and effective energy E e f / , mov- 
ing in a one-dimensional potential V e ff(r). Since 
r should be real and positive, the physically ac- 
ceptable regions are given by those r for which 
Eeff > Veff(r), owing to the square on the left 
hand side of JTOJ). Substituting for A(r) in (fig]) , 
we arrive at the equation 



AL 2 ^ 4eM 2 L 2 



(i 



4M 2 L 2 



2„4 



L z q 



(22) 

Jrt'-t IJ*^±. y*VJ 

The shape of the trajectories can be determined 
by using (|16p to express r in the following form: 



dr 
dX 



dr dip 
dip dX 



r 2 dip 



(23) 



The two special cases of particle trajectories, 
which are important from the physical point of 
view are the radial motion and the circular motion 
of the particles. Equation dl~9l) is used to study 
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radial free fall and the stability of radial trajecto- 
ries. However, analysis of the circular motion of 
particles is more con veniently done if we i ntroduce 
the variable change (jHobson et al. 1120061 ) . 



(24) 



so that 



x dXj 



-L 2 q A u w + {AM 2 L 2 



- eg > 8 



(4eM 2 - L 2 )u 6 



E z -e 



AL 1 
6 



6 

(25) 

Using equations (|2U|) to (|2"5j) . the equation 
defining the circular trajectories of test particles 
in the field of a non-rotating charged black hole in 
a five-dimensional spacetime, is obtained as 



du\ , (i 



4M*-Xj ? / + (___ 1 



E 2 ~e A\ Ae 



L 2 6 J 6u 2 L 2 
4.1. Radial motion 



P(u). (26) 



For radial motion, ip = constant and hence L = 
0. Therefore we consider equation (|18[) in the form 



£ 2 + A(r)e, 



(27) 



to analyze the radial trajectories. 

4- 1.1. Motion of massive particles 
In case of massive particles, (|27|) gives 

r 



dr 



fx) =* 2 +Mr)=E 2 - v/ 



(28) 

On differentiating (|2"5f with respect to A and di- 
viding by 2r, we get 



4M 2 2q 4 A 



dX 2 



^ + (29) 



Choosing the affine parameter A to be the 
proper time r along the path, we find that the 



force per unit mass for massive particles will be 
attractive when 



4M 2 2q 4 
— =- > -V 



Ar 
~6 ' 



(30) 



which is the condition necessary for the existence 
of bound states. As the particle moves in the 
field of the black hole, it gains kinetic energy from 
the energy of gravitational interaction. The cor- 
responding change in the gravitational potential 
energy of the particle can be estimated by consid- 
ering that the particle starts from rest at a position 
where the radial coordinate is R. Writing equation 
(|28p in terms of the proper time r of the particle, 
and assuming that at r = _R, dr/dr — 0, we ob- 
tain the kinetic energy per unit mass gained by the 
particle in terms of the change in its gravitational 
potential energy, as follows: 



i 



i 

fi 2 



2 



1 



1 

i? 1 



A 



+ §<r 2 -* 2 )- 



(31) 



Fig. [2] shows the plot of V e // vs the specific ra- 
dius r/M of the orbits, for different values of the 
specific charge q/M of the black hole, where we 
have 

for the radial motion of massive particles. 

The variation of the effective potential for 
different values of A with q/M = l/v2 is 
shown in Fig. [3J The variation is similar to the 
Schwarzschild illustrated in (IHackmann et al. 

2008b|). The similarity is due to the fact that 
the ReissnerNordstrom metric differs from the 
Schwarzschild metric only in the definition of the 
horizon function A (on account of the charge of 
the black hole) and the degree of the polynomial 

nO- 
Usrng (HU) and (HD we find that 




E l - 1 



2M 
r 



1/2 



~6~ 



(32) 
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Hence, at large distances from the black hole, 
the above derivative will blow up due to the con- 
tribution of the A-term. Thus the particles will 
not be able to reach there under the effect of the 
black hole field. However, at finite distance from 
the source, both the time derivative and the en- 
ergy of the particles must be finite, with the exact 
trajectories on the (t, r)-plane being determined 
by this equation. 

4-1-2. Motion of photons 

The equation for the radial motion of photons 
is given by 



where E is defined by (Ti4"]) . Writing 



(33) 



dr 
~dl. 



we have 



dr 



J t )=±A(r) = T 



dr dX 
d\~dt 



Ar' 

IT 



(34) 

The trajectory of the particles can be represented 
on the (t, r)-plane after integrating the above 
equation. The solution involves inverse functions. 
However, to determine the critical points and the 
nature of the particle trajectories, we have to use 
the dynamical systems analysis. The trajectories 
will also be determined from this analysis, which 
will be taken up in the next section. 

4.2. Circular motion 

For equilibrium circular orbits, u — constant 
in ([251) and hence P(u) = 0, as well as P'(u) = 0. 
From (|2l)f we have 



'AcM 2 



- 1 it 



E 2 ~c A 



L 2 6 / 6m L 



Ac 



(35) 



Thus 



P'{u) = ~6q 4 u 5 + A (^4M 2 -jjlu 3 

„ [AcM 2 \ 2 Ac 

2(^^ 1 1 u- . „,„ ■ (36) 



V L 2 



From (f36|) . applying the condition P'[u) = 0, 
we obtain the expression for the angular momen- 
tum L of a particle moving in a circular orbit of 
radius r = 1/u in RN-AdSs as follows: 



L 2 



c[6q 4 u 6 - 12M 2 



■A/2] 



3u 4 [-3<A 4 + 8M 2 u 2 - 1]' 



(37) 



Further, the condition P(u) = for the occur- 
rence of circular orbits, gives us the energy E of 
a particle executing circular motion from equation 
(j3%f as 

2 = e[2/( M )-A][/( M )-A] 

18u 4 [3g 4 M 4 -8M 2 w 2 + l]' 1 ' 

where f(u) = 3q 4 u 6 - 12M 2 u 4 + 3u 2 . 

4-2.1. Motion of massive particles 
For massive particles, (1571) becomes 



L 2 



[6g 4 u 6 - 12M 2 u 4 + A/2] 
3w 4 [-3g 4 w 4 + 8M 2 u 2 - 1] 



(39) 



Similarly, we can obtain the energy of massive par- 
ticles moving in circular orbits as 



E 2 



[2/( u )-A][/( M )-A] 
18u 4 [3g 4 u 4 - 8A/ 2 u 2 + f ] ' 



(40) 



To calculate the radius of bound circular orbits, we 
need to solve a polynomial of degree 8 in u, which 
yields eight possible values of r/M = (1/uM). For 
an orbit of a given radius, the energy E is a con- 
stant of motion. However, E varies as the radius 
of the orbits vary. The variation of E 2 and L 2 /M 2 
as a function of r/M for circular motion of massive 
particles on the equatorial plane of the black hole, 
is shown in Fig. [4] for different values of q/M. 

The graphs indicate that the energy decreases 
as the value of r/M increases for a given value of 
q/M. For a black hole with lesser value of spe- 
cific charge, particles moving in orbits closer to 
the black hole, will have higher energy and angular 
momentum. However, for sufficiently large values 
of r/M, the energy and angular momentum tends 
to a constant value. For the case of 4-dimension, 
we recall that the circular orbits remain bound 
over a range of values of r/M up to infinity, with 
the total energy tending to saturate to a stabl e 
magnitude as r/M — ¥ oo ( Hobson et al. 2006). 



6u 3 L 2 ■ 



Thus the behavior is similar even in 5-dimensions. 
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J^.2.2. Motion of photons 

The energy and angular momentum of the pho- 
tons cannot be determined from equations (|3"T)l 
and (155)) , since e = for photons. Substituting 
e = in (l26l) we obtain 



(41) 

Therefore, 

Q'(tt) = -6q 4 u 5 + 16M 2 u 3 - 2u. (42) 
For Q'(u) — 0, we have 



4M 2 ± v/!6M 4 - Sq 4 



3<Z 4 

and from Q(u) — we obtain 



(43) 



|^ = q 4 ,i 6 - 4M V + u 2 - £ . (44) 

6 



Equation (|43|) indicates that the photons can have 
circular trajectories for only two specific values of 
r/M for a black hole of a given mass and charge. 
This feature is an important distinction between 
the circular motion of photons and those of mas- 
sive particles discussed above. 

4.3. Stability of orbits of massive particles 

The stable circular orbits occur for those values 
of r which are located at the local minimum of 
the potential. The local maxima of these curves 
correspond to the radii of unstable circular orbits. 
Fig. [5] shows the plot of the effective potential vs 
r/M for different values of q/M. 

Fig. IS] shows the variation of the effective po- 
tential for different values of A with q/M = 1.1. 

The minimum radius for stable circular orbits 
will occur at the point of inflexion of the function 
P(u), for which we must have P"(u) = 0, along 
with P(u) = and P'{u) = 0. From ((36j), we get 



P"(u) 



4„,4 



-30<fu 



48 AT 



8M 2 



- 2 



12 g 4 
L 2 



A 



IP- ~J>W ^ 

Thus to find the radius of the innermost stable 
circular orbit of massive particles, we have to solve 



a polynomial of degree eight in u. However, we can 
reduce it to a polynomial of degree 4 with the help 
of the substitution u = *Jx. This leads us to the 
equation 

30L 2 q 4 x 4 -{4:8M 2 L 2 -l2q 4 )x 3 -(8M 2 ~2L 2 )x 2 -A = 

(46) 

The radius of the innermost stable circular orbit 
of massive particles in a five-dimensional RN-AdSs 
can be calculated from (|46l) . provided 4M 2 L 2 =^ q 4 
and 4M 2 =t L 2 . These tracks will be visible to an 
observer only if this minimum radius is greater 
than the radius of the event horizon. 



4-3.1. Special case 

For 4M 2 L 2 = q 4 , 4M 2 ^ L 2 since A < 0, we 
find that 



r ISCOM 



3QLq z 



(4M 2 - L 2 ) + ^/(4M 2 -i 2 ) 2 +30LVA' 

(47) 

indicating that in this case there will be only two 
possible real values of the radius, provided (4M 2 — 
i 2 ) 2 > |30LVA|. 

4.4. Stability of orbits of photons 

To analyze the stability of photon trajectories, 
we consider the first integral of the geodesic equa- 
tions 

fdr\ 2 ., ..L 2 



which can be recasted as 

2 



1 fair 
I 2 \ dX 



+ V ef f(r) = E eff , 



(48) 



(49) 



with 



V, 



eff 



and 



4M 2 



E 

Eeff - 



Ar 2 



Making the transformation A — > LX we can 
rewrite this equation in the form of an equation 
for effective potential V e /f. Fig. [7] shows the plot 
of V ef f vs r/M. 
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When the photons describe circular trajecto- 
ries, the minimum radius r mc of the stable cir- 
cular orbits can be obtained from the condition 
Q"(u) = 0, where 

Q"(u) = -30 9 4 u 4 + 48M 2 m 2 - 2. (50) 

Thus we get 



15cT 



12M 2 - \/l44Af 4 - 15g 4 



5. Analysis of particle trajectories using 
the dynamical systems approach 

In this section, we consider the geodesic equa- 
tions as a dynamical system and proceed to deter- 
mine its phase space. A number of authors have 
used this method in various context ( Uzan et al 



20011 iDahia et al. |[2007tlWainwright et al. Il997j l 

For the analysis of the phase trajectories and 
the nature of cr itical points, we u se the method 
adopted earlier ( Guha et al. Il2010h . 

Let us define the dynamical system in terms of 
three new variables as follows: 



V=- 
dA' 



and 



W 



chp_ 
dX' 



(51) 
(52) 

(53) 



The geodesies equations (1101) to ([12]) , when written 
in terms of these variables, reduce to the following 
form: 



^ + A{r)B{r)U 2 - ^-V 2 + rA{r)W 2 = 0, 
dX 



and 



A(r) 
dW 1. 



-VW = 0. 



(55) 
(56) 



d\ r 

These three variables are related through the first 
integral of geodesic equations i.e. 

A(r)U 2 - -^-V 2 + r 2 W 2 = -e. (57) 
v ; A{r) 



It is not necessary to study the complete system 
(|54p to (|56|) . since W can be determined alge- 
braically from ([5T]). Substituting for W from (|57|) 
in ((55)) . we obtain 



f + f(rfl(r)-A(rr + ^(r) 



rB{r))V 2 



eA{r) 



0. (58) 



The description of the phase trajectories now re- 
duces to the analysis of the real non-linear dynam- 
ical system defined by the equations 



where 



and 



± = V 
dX 

dU 

fx = J{ ' J ' v ' r) ' 



B(r) 

H(U,V,r) = --^-UV 
A(r) 



(59) 



(60) 



J(U, V, r) = - ^l(rB(r) - A(r))U 2 



r r 
1 



rA(r) 



(A(r) - rB(r))V 2 



with e = 1 , for timelike and null geodesies re- 
spectively. Both H and J have continuous first 
partial derivatives for all (U,V). 

5.1. Determination of the fixed points 

The fixed points of the above dynamical system 
are the solutions of the system (|59l) , for which V — 
0, H = and J = 0. These solutions represent the 
equilibrium positions for the particle trajectories, 
which can be either stable or unstable. 

Let {Uq.Vq = 0, 7"o) represent the fixed point of 
the phase trajectories. We can see that Vo = is 
sufficient to satisfy H = 0. It also implies that r 
is constant for the fixed point. To have J = in 
this condition, we further need A(r) = —f(r) = 0, 
for both timelike and null geodesies, where f(r) 
is given by equation ((4}. In fact, the value of r 
can be determined from the condition A(r) = 0, 
which means that tq is constant for a black hole 
of a given mass M and charge q, when A remains 



S 



constant. Consequently, for a given black hole, 
the fixed points will lie on a definite phase plane 
corresponding to a constant value of ro. Finally, 
the coordinate Uq can be evaluated from 



eA(r ) A(r ) 



ro 



(r S(r ) - A{r ))U 2 = 0, (61) 



where e = 0, 1 respectively for null geodesies and 
timelike geodesies. Thus we have 



(62) 



(r B(r ) - A(r Q ))- 

For null geodesies e = 0, so that 

U = 0. (63) 

Hence, for a black hole with a given M and q, (0,0) 
is the fixed point on the (U-V) phase plane defined 
by ro = a where a is a constant. 
For timelike geodesies, e = 1 and 



U a = ±< 



1 



(r B(r ) - A(r )) ' 



(64) 



Evidently, the timelike geodesies will possess defi- 
nite fixed points for a black hole of a given M and 

q- 

5.2. Analysis of phase trajectories 

The phase evolution of the system in the (U, V) 
phase plane is given by the solution of the set (|59p . 
We note that the equation 



dV 
dU 



J_ 

H 



(65) 



specifies the phase evolution of the system in the 
(U, V) phase plane, provided H ^ at this point, 
which is true except at the critical point (0,0) for 
the null geodesies. Substituting the expressions 
for H and J in ([65|) and integrating, we obtain 



A(r)C(r)U 3 eA(r)U C{r)UV 2 B(r)UV 2 



3r 



rA(r) 



2A(r) 



D = 0, 



(66) 

where D is a constant of integration and C(r) = 
rB(r) — A(r). We can rescale our coordinates so 
as to reduce this constant D to zero and simplify 
to arrive at the result 



fiU - hV 2 - h = 0, 



(67) 



where /i, f% and are functions of r. For null 
geodesies, we have /a = 0. 

Below Fig. E] and Fig. [9] show possible trajec- 
tories of massive particles for different choices of 
parameters /i , /2 and f% . We find that the trajec- 
tories are either elliptic or a pair of hyperbolas. 

The phase trajectories of the photons are a pair 
of straight lines intersecting at (0,0), with their 
slopes changing according to the choice of the pa- 
rameters /i and /2. Fig. [TOl shows a sample plot 
of the null geodesies. 

5.3. Nature of orbits 

Our investigations using the method of dynami- 
cal systems, indicate that the geodesies of massive 
particles near a non-rotating charged black hole 
in a five-dimensional spacetime possesses definite 
fixed points for a given black hole and are either 
periodic bound or escape orbits. Moreover, for a 
black hole of a given charge and mass, the null 
geodesies possesses a unique fixed point (U = 0, 
Vq = 0, ro), which correspond to the condition 
4r = and — along the geodesic curves. 
Hence the null geodesies are terminating orbits. 

6. Summary and conclusions 

We have studied the timelike and null geodesies 
of neutral particles in the field of a non-rotating, 
charged black hole in a five-dimensional spacetime 
from the point of view of effective potential for- 
malism and the dynamical systems approach. The 
black hole horizons are located with the help of the 
plot of the effective potential. In addition to ra- 
dial trajectories, we have investigated the circular 
trajectories for photons and massive particles, as 
well as the fixed points of particle trajectories. It is 
found that photons will trace out circular trajecto- 
ries for only two distinct values of specific radius 
of the orbits. The radius of the innermost sta- 
ble circular orbit of massive particles is totally de- 
fined in terms of their angular momentum and the 
mass and charge of the black hole. To determine 
the nature of tra jectories and the fixed points, we 
have used the dynamical systems analysis, defining 
the dynamical variables in terms of the derivatives 
of the coordinates with respect to the affine pa- 
rameter along the geodesic curves. Consequently, 
we have found that the geodesies of massive par- 
ticles near a non-rotating charged black hole in 
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5D are either periodic bound or escape orbits and 
there are definite fixed point for these trajectories. 
Moreover, the null geodesies have a unique fixed 
point (Uq — 0, Vq — 0, 7*0,) and these orbits are 
terminating orbits. 

We conclude with the note that here we as- 
sumed the singularity to be a black hole and not a 
naked singularity. We know that black holes and 
naked singularities can be observationally differen- 
tiated through their gravitational lensing features. 
Whether these singularities can be differentiated 
in terms of the nature of the particle trajectories 
is the subject of our investigation and will be re- 
ported in a future work. 
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Fig. 1. — Diagram showing the plot of effective 
potential V e ff vs radial distance r/M in the case 
of radial motion of massive particles with q/M = 
0.5 for different range of values of r/M for sample 
values of A = -0.05 (red) and A = 0.05 (blue) in 
the top figure and A = —0.0005 (red) and A = 
0.005 (blue) in the bottom figure 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 3. — Diagram showing the plot of effective 
potential V e ff over different range of values of 
r/M for sample values of A = —0.5 (orange) 
A = -0.0005 (black) and A = 0.1 (green) in the 
case of radial motion of massive particles having 
q/M = 1/V2 
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Fig. 4. — Diagram showing the radial dependence 
of E 2 and L 2 /M 2 for particles moving in circular 
orbits on the equatorial plane with q/M = 0.5 
(black line) and q/M = 1.414 (blue line) 



-1.5 



Fig. 5. — Diagram showing the variation of the ef- 
fective potential V e f / with a variation of the spe- 
cific radius of circular trajectories of massive par- 
ticles for increasing values of q/M with q/M =1.1 
(black), q/M = 1.23 (blue) and q/M = 1.414 (red) 
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Fig. 6. — Diagram showing the plot of effective 
potential V e ff over different range of values of 
r/M for sample values of A = —0.05 (orange), 
A = 0. - 0005 (black) and A = 0.005 (green) in 
the case of circular trajectories of massive particles 
having q/M = 1.1 



Fig. 7. — Variation of the effective potential V e ff 
with the variation of the specific radius of circular 
trajectories after the transformation from A to LX. 
On the left, A = -0.0005 with q/M = 1 (blue) and 
q/M = 1.19 (green). On the right, q/M = 1 with 
A = -1.0 (orange) and A = +1.0 (black) 
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